Comparative putative metabolites profiling of Tachypleus gigas and Carcinoscorpius rotundicauda hemocytes stimulated with lipopolysaccharide

Horseshoe crabs are among the most studied invertebrates due to their unique, innate immune system and biological processes. The metabolomics study was conducted on lipopolysaccharide (LPS)-stimulated and non-stimulated hemocytes isolated from the Malaysian Tachypleus gigas and Carcinoscorpius rotundicauda. LC–TOF–MS, multivariate analyses, principal component analysis (PCA), and partial least squares-discriminant analysis (PLS-DA) were included in this study to profile the metabolites. A total of 37 metabolites were identified to be differentially abundant and were selected based on VIP > 1. However, of the 37 putative metabolites, only 23 were found to be significant with ANOVA at p < 0.05. The metabolites were identified using several databases, and the literature review of the metabolites was reported in the manuscript. Thus, this study has provided further insights into the putative metabolites' presence in the hemocytes of horseshoe crabs that are stimulated and non-stimulated with LPS and their abundance in each species. Several putative metabolites showed they have medicinal values from previous studies.


Extraction of metabolite
The extraction was optimized based on Muelas et al. 15 .First, approximately 200 μl of cold extraction solvent [methanol/acetonitrile 50/50 (% v/v)] was added to 600 μl hemocytes stimulated and non-stimulated samples.The samples were mixed vigorously by vortex for approximately 30 s.After vigorously mixed, the samples were incubated at − 20 °C for 20 min.The samples were centrifuged at 13,000 rpm for 5 min at 4 °C.After being centrifuged, the supernatant was transferred and dried using a speed vacuum.The methanol extracts were stored at − 80 °C until the sample analysis.

Analysis of LC-TOF-MS
The mass spectrometry analysis was performed at the Institute of Systems Biology (INBIOSIS), Universiti Kebangsaan Malaysia.Three biological and six technical replicates from each group were analyzed using LC-TOF-MS.The analyses were performed using ultra-high-performance liquid chromatography (UHPLC) with a microTOF Q III mass spectrometer (MS) (Bruker Daltonics, Bremen, Germany) equipped with an electrospray source (ESI) and connected to an Ultimate 3000 UHPLC system (Dionex, Sunnyvale, CA, USA) equipped with an Acclaims Polar Advantage II, 3 × 150 mm, 3-μl particle size C18, reverse-phase column at the flow rate of 0.4 ml/min at 40 °C with a sample injection volume of 3 μl.The running buffers were (A) deionized distilled water with 0.1% formic acid and (B) 100% acetonitrile with a 22-min total run time.The gradient elution was carried out from 5 to 80% buffer B over 22 min as follows: 5% B (0-3 min), 80% B (3-10 min), 80% B (10-15 min), 5% B (15-22 min).MS was performed in ESI positive ionization mode with the capillary voltage at 4,500 V, nebulizer pressure at 1.2 bar, and drying gas at 8 L/min in 200 °C.The scan range was from 100 to 1000 m/z.Data processing was performed using the Data Analysis 4.0 and Profile Analysis (Bruker Daltonics) software.Lastly, the graphical representation was performed on WEGO for data analysis of the sample.The caffeic acid standard was used as a control to ensure data quality by providing the optimum condition of extraction and instrument.

Statistical analysis
Multivariate statistical evaluation of the preprocessed metabolic profiling data was performed with SIMCA-P + (version 12) (Umetrics, Umea, Sweden).Bucketing was generated using Data Analysis 4.0 and Profile Analysis (Bruker Daltonics).The differential metabolites obtained from the multivariate data analysis software were validated using ANOVA with post hoc Tukey's tests and t-tests.All metabolites contributing to group separation were significant at p < 0.05.The Permutation test was performed for validation.

Compound identification
Compound identification of metabolites was performed by comparing the accuracy of the m/z value < 10 ppm and MS/MS spectra with available online databases: Human Metabolome Database (HMDB) 16 , Metabolite and Chemical Entity Database (METLIN) 17 , KEGG 18 and LIPIDMAPS 19 and MassBank.KEGG pathway analysis was used to determine the putative metabolite pathway.
www.nature.com/scientificreports/O-antigen, and hydrophilic core polysaccharide 20 .It is a primary factor in hemocyte activation and is regarded as one of the most important experimental tools for understanding the horseshoe crab immune response.
The metabolomics studies used a high throughput method (LC-TOF-MS) and multivariate analyses (PCA and PLS-DA).Based on the differences in the peak and area of the chromatogram observed following LC-TOF-MS analysis in Fig. 1, there are differences in the putative metabolite profile in both species stimulated and non-stimulated hemocytes.Multivariate analyses were then conducted to get a general overview and understanding of the spread of variability in the data.
This study used PCA and PLS-DA to analyze the preprocessed LC-TOF-MS datasets.The PCA model reveals the general metabolic information and visually eliminates abnormal sample data.It was also conducted to determine the global differences between the metabolic profiles of the groups.Based on the analyses, all samples from the two species appeared in the Hotelling T2 with 95% confidence, suggesting that all the samples can be used for further research.The parameters described the PCA model (R 2 X = 0.0918, Q 2 = 0.246).
As depicted in Fig. 2, the PCA plot shows four different groups of the samples, which are C. rotundicauda stimulated with LPS (green), non-stimulated C. rotundicauda (blue), T. gigas stimulated with LPS (red) and non -stimulated T. gigas (yellow).Based on the plot, there are differences between both horseshoe crab species, which can be measured from the locations of the metabolite's distributions in the quadrant.The C. rotundicauda (green and blue) samples were clustered on the lower quadrant on the right side, while the T. gigas (yellow and red) samples were clustered on the left upper and lower quadrants.
Comparison between LPS-stimulated (red) and non-stimulated hemocytes (yellow) of T. gigas, reveals an apparent difference in the metabolites produced.In contrast, in C. rotundicauda, there is no evident difference in the metabolites produced by the hemocytes stimulated with LPS (green) and without stimulation with LPS (blue), as they were grouped within the same quadrant.
The PLS-DA model analysis further demonstrated distinct discrimination in the metabolomic changes between the two species (Fig. 3).The acceptable values for the intercepts R 2 (cum), the goodness of fit, was 0.99, and Q 2 (cum), predictability, was 0.78.Model cross-validation through permutation tests (100 permutations) and sevenfold cross-validation generated the intercept R 2 and Q 2 0.978 and 0.267, respectively.The results show that the PLS-DA model is not overfitting and is valid for this metabolomic profiling.
The data presented herein demonstrated a clear and significant separation between the two species of horseshoe crabs and significant differences in the putative metabolites of the stimulated and non-stimulated hemocytes of T. gigas by multivariate analyses using PCA and PLS-DA.
PLS-DA also allows for determining and discriminating metabolites using the variable importance on projection, known as VIP.The VIP score value indicates the contribution of a variable to the discrimination between all the classes of samples.Mathematically, these scores are calculated for each variable as a weighted sum of squares of PLS weights.The mean VIP value is 1, and VIP values over one are usually considered significant.A high score agrees with a solid discriminatory ability and thus constitutes a criterion for selecting biomarkers.The discriminating metabolites were obtained using a statistically significant threshold of VIP values obtained from the PLS-DA model on the normalized raw data at the univariate analysis level.The P value was calculated www.nature.com/scientificreports/by one-way analysis of variance (ANOVA) for four groups analysis.Metabolites with VIP values greater than 1.0 and p-values less than 0.05 were considered statistically significant metabolites.Figure 4 depicts a PLS-DA loading plot showing the variables that contributed to separating the four groups.
Based on Fig. 4, the farther the point is from the origin, the greater the weight value or the greater the effect of determining the grouping of the samples.In the non-stimulated C. rotundicauda hemocytes, the metabolite masses that distinguished it from the other three groups are 6.

Compound identification
Identification of the compounds was performed by comparing the accuracy of the m/z value (< 10 ppm) and MS/ MS spectra with available online databases: Human Metabolome Database (HMDB) 16 , Metabolite and Chemical Entity Database (METLIN) 17 , KEGG 18 , LIPIDMAPS 19 and MassBank.All the putatively identified metabolites are reported in Table 1.

Comparative metabolomics analyses
An untargeted metabolite profiling of two horseshoe crab species, C. rotundicauda, and T. gigas, was performed in this study.There were two conditions observed: stimulated and non-stimulated with LPS.Multivariate statistical analyses such as PCA, PLS-DA, ANOVA, Post Hoc Tukey test and t-test were conducted.Based on the PLS-DA analysis, 37 metabolites are in the VIP > 1 category.Of 37 metabolites, only 23 were statistically significant when ANOVA was performed at p-value < 0.05.All the putative metabolites were tabulated in Table 1.Further statistical analysis using Post Hoc Tukey was also conducted after performing ANOVA.Of 23, 7, 8-Dihydroxanthopterin, Fig. 5a is one of the putative metabolites found to be significant in both tests.No research studies have reported on this metabolite in horseshoe crabs.However, several studies on this metabolite from other organisms have been found.For example, in Stizostedion lucioperca, this metabolite influences the eye colour and vision of the fish 21 .The guanine crystals form a reflective layer that produces the silvery colour present on the eye surface.Meanwhile, the block-shaped crystals backscatter light into the retina, which helps to increase the sensitivity to light.
In human studies, this putative metabolite was found in the urine of phenylketonuria and lethal hyperphenylalaninemia patients 22 .It was reported to play a role in the pathogenesis of neurological symptoms in both diseases.Figure 5a shows that this putative metabolite was reduced in both species' hemocytes after LPS stimulation.However, Tukey's multiple comparisons test shows that the reduction of metabolites after stimulation with LPS is not statistically significant in C. rotundicauda but statistically significant in T. gigas.A comparison between species using Tukey's also shows significant differences between both species, and this metabolite can be found more abundantly in C. rotundicauda than in T. gigas.
Next, Harman, Fig. 5b, a natural B-carboline alkaloid becoming interesting due to its anti-cancer properties 23 .It also was found to be decreased in C. rotundicauda and T. gigas after hemocytes were stimulated with LPS.However, despite reductions, Tukey's multiple comparisons test shows the reduction is insignificant in each species between stimulated and non-stimulated forms.Despite that, there are significant differences in the metabolite between both species in which this metabolite is abundant in C. rotundicauda.By knowing this, isolation of the metabolites can be suggested on C. rotundicauda instead of T. gigas.This metabolite is usually known to be derived from plants, and sesame seed oil was reported to have high levels of β-carbolines 24 .The fungi entomopathogen Conidiobolus coronatus 25 can also secrete it.It is also considered a nonpolar heterocyclic aromatic amine with potential mutagenicity 24 .It is also a reversible competitive monoamine oxidase inhibitor, increasing serum serotonin concentrations in tissues 25 .Antibacterial activities of Harman analogues against four gram-positive and two gram-negative bacteria damaged bacterial cell membranes and walls and disrupted the function of type II topoisomerase 26 .These derivatives also have potential as new bactericides and antibiotics, as the in-vivo antibacterial assay shows a protective efficacy of 81% 26 .In insects, Harman also resulted in delayed pupation and adult eclosion and inhibited total monoamine oxidase activity 25 .
Phosphatidylglycerol, Fig. 5c, increased in hemocytes after stimulation with LPS in both species.However, despite having an increment, Tukey's multiple comparisons test shows the increment is insignificant in both species stimulated and non-stimulated groups.Despite that, there is a significant difference if comparing two species in which this metabolite can be found to be abundant in T. gigas.Phosphatidylglycerol (PG) is a naturally occurring phospholipid and is essential for the growth and photosynthesis of photosynthetic organisms 27 .It is the only major phospholipid in the thylakoid membrane of chloroplasts 27,28 .As it is crucial for photosynthesis, Vol:.( 1234567890 www.nature.com/scientificreports/ the loss of PG in Arabidopsis thaliana resulted in severe defects in the growth and development of chloroplast with decreased accumulation of chlorophyll, impaired thylakoid formation, and also downregulation of photosynthesis-associated genes encoded in nuclear and plastid genomes 28 .PG is also one of the components needed in daptomycin to exert its antibacterial effect 29 .PG and sulfoquinovosyldiacylglycerol (SQDG) have similar physicochemical properties, bilayer thickness, and bending rigidity 30 .However, the function of this metabolite in horseshoe crabs should be further elucidated.Androstenedione, Fig. 5d increased after the stimulation of hemocytes with LPS in both species.However, the Tukey test shows that this increment is not statistically significant.Despite that, both species have a significant difference in metabolites.In C. rotundicauda hemocyte, this metabolite is abundant compared to in T. gigas.Thus, if Androstenedione isolation is yet to be performed, it can be isolated from C. rotundicauda.Several studies have reported on androstenedione, but no specific analysis of androstenedione of horseshoe crabs has previously been reported.Dong-Ma et al. reported that androstenedione and androgens androstenedione (ADD) are predominant steroid hormones in surface water or wastewater and can disrupt the endocrine system in fish 31 .Androstenedione is produced in male and female gonads and the adrenal glands and is known for its crucial role in producing estrogen and testosterone 32 .It is also a precursor for several steroid substances like testosterone, estradiol, ethinyl estradiol, testolactone, progesterone, cortisone, cortisol, prednisone, and prednisolone 33 .It is also sold as an oral supplement to increase testosterone levels 32 .The supplement can also lower triglycerides (TG) and high-density lipoprotein (HDL) cholesterol, increase oestradiol concentration, and is a natural alternative to an anabolic steroid 32,34 .Other uses of the metabolite include as an enhancer for athletic performance, building body muscles, reducing fats, increasing energy, maintaining healthy RBCs, and increasing sexual performance 32 .Androstenedione is also listed among performance-enhancing drugs (PEDs).However, it was banned by the World Anti-Doping Agency and International Olympic Committee 32 .
Diglycerides, Fig. 5e, were found to be increased after stimulation with LPS in C. rotundicauda; meanwhile, in T. Gigas, the stimulation of hemocytes with LPS decreased the metabolites abundantly.Statistical analysis using Tukey shows that the changes are significant in T. gigas but not C. rotundicauda.In a study conducted by Song et al. 35 , monoglyceride and diglyceride were shown to have antiviral and antibacterial properties and act as emulsifiers to increase the digestibility of dietary lipids.Its supplementation could also effectively reduce fat loss, decrease inflammatory factor levels, and control total cholesterol concentrations during lactation 35 .In brown adipose tissue, L-Carnitine helps increase TG and diglyceride levels and reduces glycerophospholipids and sphingolipids.
L-Iditol, Fig. 5f increased after stimulation with LPS in both species.Statistical analysis using Tukey shows that only changes in T. gigas are significant.Studies regarding L-iditol alone are underreported.Only two relevant studies have been reported regarding L-Iditol, the angiosperm-Yunnanopilia longistaminata, a new plant source for L-iditol and taxanes 36 .Secondly, a series of quaternary diammonium salt derivatives of 1,4:3,6-dianhydrol-iditol were synthesized, and two quaternary ammonium salts (QAS) with octyl and decyl residues exhibited antimicrobial activity 37 .
Octaprenyl diphosphate Fig. 5g also was found to be decreased after stimulation with LPS in hemocytes of both species and is statistically significant in the Tukey test.This metabolite is essential for the normal growth of Escherichia coli 38 .However, the relation between this metabolite and horseshoe crab is unknown.
Valine, Fig. 5h, was found to be increased in C. rotundicauda; meanwhile, in T. gigas, it was found to be decreased after stimulation with LPS.However, the changes were not statistically significant when a Tukey comparison was conducted.Despite that, this metabolite is abundant in T. gigas, which means it can be isolated in this species if needed.Valine has extensive industrial applications and is an intermediate for synthesizing agricultural pesticides and semisynthetic veterinary antibiotics 39 .Bacillus cereus could have a potential for industrial www.nature.com/scientificreports/production of valine under optimized conditions 39 .Dietary L-Valine supplementation modulates the inflammatory response and microbial metabolites 40 .Anthranilate Fig. 5i increased in C. rotundicauda and decreased in T. gigas after hemocyte stimulation with LPS.However, analysis using Tukey shows that the metabolite changes are not significant.Despite that, this metabolite was found to be abundant in T. gigas.This metabolite is widely used as a precursor in producing dyes, fragrances, plastics, and pharmaceutical compounds 41,42 .Microorganisms produce Anthranilate as an intermediate in the tryptophan biosynthetic pathway 41 .It has various biological activities, such as anti-inflammatory, antineoplastics, anti-malarial, and has α-glucosidase inhibitory properties 43 .Methyl anthranilate (2-aminobenzoic acid methyl ester) irritates birds' senses of taste and smell 44 , protecting sweet cherry orchards against birds.Anthranilate also increased the antibiotic susceptibility of other species of bacteria, such as Escherichia coli, Salmonella enterica, Bacillus subtilis, and Staphylococcus aureus 45 .Evaluating the antifungal activity in vitro of the active films containing methyl anthranilate showed great effectiveness against Penicillium expansum and Botrytis cinerea, demonstrating the potential applicability of the developed films for active food packaging.Evodileptin B (1) is a natural anthranilate derivative isolated from the ethanol extract of the aerial parts of Evodia lepta (Spreng.)Merr., a traditional medicinal plant of the family Rutaceae 46 .Evodileptin B has solid neuroprotective properties and may help treat Parkinson's Disease 46 .Linalyl anthranilate (LNA) generates reactive oxygen species, initiates lipid peroxidation, and damages the bacterial membrane, resulting in intracellular leakage and eventually killing Klebsiella pneumoniae 47 .Another study of Anthranilate is a novel anthranilate analogue (SI-W052) that inhibited LPS-induced tumour necrosis factor (TNF)-α and interleukin (IL)-6 on microglia 48 .However, further studies must be conducted to understand its relationship with horseshoes crab metabolites.
N-acetyllactosamine Fig. 5j was found to be decreased in C. rotundicauda, after being stimulated with LPS, but it increased in T. gigas.However, the Tukey test shows that the increment and reduction are not statistically significant.N-acetyllactosamine (LacNAc), specifically β-d-galactopyranosyl-1,4-N-acetyl-d-glucosamine, is a unique acyl-amino sugar and a critical structural unit in human milk oligosaccharides, an antigen component of many glycoproteins, and an active antiviral property for the development of effective drugs against viruses 49 .The 6-sulfo -N-acetyllactosamine was found to inhibit the binding of the SARS-CoV-2 spike protein S1 subunit with blood group A RBCs and reduce the interaction between the spike protein S1 subunit and Angiotensinconverting enzyme 2 (ACE2) in SARS-CoV-2 infection 50 .
We performed the t-test to analyze each species' pre-and post-LPS treatment metabolite changes.Only three metabolites exhibited a significant increase-specifically, the increments in L-iditol and Diacylglycerophospholipid in T. gigas.Meanwhile, dtMP deoxythymidylic acid was statistically significant in C. rotundicauda (Table S1).
In summary, the significant changes observed in metabolites by the t-test, such as L-iditol, Diacylglycerophospholipid, and dtMP deoxythymidylic acid, following LPS stimulation in T. gigas and C. rotundicauda, suggest their involvement in crucial pathways; immune regulation, inflammation, and stress responses.These metabolites likely play roles in modulating immune cell function, membrane stability, or nucleotide metabolism, contributing to the immune response.The observed differences in metabolite abundance between the two species indicate unique strategies in response to the immune challenge posed by LPS.Species-specific responses, exemplified by T. gigas' significant decreases in 7,8-Dihydroxanthopterin, Phosphatidylinositol PI, and octaprenyl diphosphate, and C. rotundicauda's specific decrease in octaprenyl diphosphate, highlight distinct molecular adaptations which evolutionary history, environmental factors, and ecological niches might influence.These species-specific strategies reflect C. rotundicauda and T. gigas' adaptations to diverse environments and immune challenges.

Pathway enrichment analysis
MetaboAnalyst 51 shows there are three significant pathways where the metabolites are found to be enriched, which are glycerophospholipid metabolism, valine, leucine, and isoleucine biosynthesis and glycosylphosphatidylinositol (GPI)-anchor biosynthesis (Fig. 6).Glycerophospholipids are the most abundant and dominant in cell membranes as they provide stability, fluidity, and permeability 52 .Moreover, they must function correctly as www.nature.com/scientificreports/membrane proteins, receptors, and ion channels and as reservoirs for second messengers and their precursors.Thus, phosphatidylglycerol (12.46 min: 779.5176 m/z) in horseshoe crabs probably helps with cell membranes' excellent structure and signalling.On the other hand, the biosynthesis of valine, leucine, and isoleucine is crucial as they play critical roles in regulating energy homeostasis, nutrition metabolism, gut health, immunity, and disease in humans and animals 53 .Perhaps the existence of the metabolite valine (1.82 min: 140.0685 m/z) in horseshoe crabs also plays the same role as in other organisms.Glycosylphosphatidylinositol functions as an anchor to link cell membranes and proteins.These proteins act as enzymes, adhesion molecules, complement regulators, or co-receptors in signal transduction pathways 54 .The richness of this metabolite in this pathway is probably linked to the first hit, the glycerophospholipid metabolism, as it seems to work together.However, the metabolite name (1-Phosphatidyl-D-myo-inositol) that hit this pathway is not precisely the same as our putative metabolite Phosphatidylinositol (12.26 min: 911.5912 m/z).Despite that, as we refer to the KEGG database, it refers to the same metabolite but differs in name.The KEGG pathways involved are map00564 Glycerophospholipid metabolism, map00290 Valine, leucine and isoleucine biosynthesis and map00563 Glycosylphosphatidylinositol (GPI)-anchor biosynthesis [55][56][57] .

Transcriptomics, proteomics, and metabolomics of hemocytes after LPS stimulation
Previous studies had conducted transcriptomics and proteomisc using hemocytes of horseshoe crabs 11,14 .The transcriptomics analysis reported 1338 genes were significantly upregulated, and 215 genes were downregulated after hemocytes were stimulated with LPS.Meanwhile, proteomics analysis reported 154 proteins were identified in the stimulated and non-stimulated form of hemocytes.From 154 proteins, 54 were found to be unique in hemocytes stimulated with LPS, and 25 were unique in non-stimulated form.Thirty-seven proteins are found to be shared in both conditions.Tachylectin-2, coagulogen, c-reactive proteins, histones, hemocyanin, and DNA polymerase, all of which play essential roles in the organism's innate immunity, were found to be differentially expressed in hemocytes after the LPS challenge 14 .Gene ontology enrichment analysis from both studies showed several differentially expressed genes and proteins predictively involved in several metabolic processes such as cellular metabolic process, protein metabolic process, macromolecule metabolic process, and organonitrogen compound metabolic process.Indeed, our study showed several putative metabolites such as Androstenedione, 7,8-Dihydroneopterin, and Phosphatidylglycerol involved in the metabolic pathways by KEGG pathway analysis.Metabolites such as 7,8-Dihydroxanthopterin and Harman involved in the KEGG biosynthesis pathway, which the pathway was also enriched in a study by Sarmiento et al. 11 at the gene level.
Nevertheless, it requires extensive study to understand the metabolites pathway and its functions to correlate our findings with previous findings at the gene and protein levels.Despite that, several metabolites such as Harman, L-Iditol, contignasterol, valine, Anthranilate, N-acetylactosamine, and Glycerophosphoinositol could be validate in further studies as they have immense benefits.

Figure 1 .
Figure 1.The base peak chromatogram (BPC) of non-stimulated hemocyte cells of C. rotundicauda A(i) vs hemocytes stimulated with LPS A(ii) and non -stimulated hemocyte cells of T. gigas B(i) vs stimulated with LPS B(ii) by Liquid Chromatography-Time of Flight-Mass Spectrometry (LC-TOF-MS).

Figure 2 .
Figure 2. PCA score plot of PC1 versus PC2 scores for compounds or metabolites detected in each treatment: Hemocytes of C. rotundicauda and T.gigas non-stimulated and stimulated with LPS.

Figure 3 .
Figure 3. PLS-DA score plots for compounds or metabolites detected in each treatment: Hemocytes of C. rotundicauda and T. gigas non-stimulated and stimulated with LPS (C.rotundicauda stimulated with LPS (green), C. rotundicauda non-stimulated (blue), T. gigas stimulated with LPS (red), and T. gigas non-stimulated (yellow).

Figure 4 .
Figure 4. Loading plots of PLS-DA for metabolites (masses) detected via untargeted LC-TOF-MS in different treatments.The green dots labelled with retention time represent the masses distributed w*c(1) and w*c(2) planes and the highest VIP (Variables Importance for the Projection) > 1. Blue dots represent the origin.

Figure 6 .
Figure 6.Overview of enrichment significant putative metabolites of horseshoe crabs.

Table 1 .
Putatively identified metabolites with the highest variable importance on projection (VIP) score as determined by a partial least square-discriminant analysis (PLS-DA) and one-way ANOVA.